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Since late 2007, we have been monitoring a large sample of known and likely gamma-ray— loud blazars at 15 GHz 
twice per week with the Owens Valley Radio Observatory (OVRO) 40 m Telescope. Our initial sample included 
the 1158 sources above declination —20° from the Candidate Gamma-Ray Blazar Survey (CGRaBS), and we 
have since added nearly 400 more sources, including all blazars associated with Fermi Large Area Telescope 
(LAT) detections in the First AGN Catalog (1LAC). Here, we describe the new sample and present results for 
2008 through early 2011. Using statistical likelihood analyses, we compare the variability amplitude for various 
sub-populations within our sample. These include comparisons of gamma-ray-loud versus gamma-ray-quiet 
objects, BL Lac objects versus flat-spectrum radio quasars, and a study of the variability amplitude trend with 
rcdshift. We also describe KuPol, the new digital Ku-band receiver being constructed for the 40 m telescope. 
This new receiver will provide total intensity and linear polarization measurements over the 12-18 GHz band, 
with 16 MHz spectral resolution. 



1. INTRODUCTION 

Blazars are the most extreme class of active galac- 
tic nuclei (AGN), probably resulting when an AGN 
is viewed along its relativistic jet axis. Although a 
variety of models have been proposed, and despite 
decades of observation, many fundamental questions 
about blazars and AGN physics remain open. Among 
them: 

• How are jets launched, accelerated, collimatcd, 
and confined? 

• Of what are the jets composed? 

• Where in the jet are the observed photons emit- 
ted? 

• What are the details of the emission mecha- 
nisms? 

Gamma-rays are our best tool for probing the most 
extreme processes within AGN, and coordinated mul- 
tiwavelength observations are crucial to the interpre- 
tation of gamma-ray data. 

1 .1 . OVRO 40 m Program 

Since late 2007, we have undertaken a continu- 
ous, fast-cadence (twice weekly per source) monitoring 
program beginning with the systematically-selected 
sample of 1158 blazars north of S = —20° from the 
Candidate Gamma-Ray Blazar Survey (CGRaBS) [2]. 
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The sample has grown to more than 1550 sources, in- 
cluding all 454 "clean" AGN associations from the 
first-year Fermi Large Area Telescope (LAT) AGN 
catalog (1LAC) pQ. Newly LAT-detected AGN are 
added to the sample as well. Results for the CGRaBS 
sample during the 2008-2009 period are presented 
and analyzed and the program is described in de- 
tail in [6]. Radio light curves for the CGRaBS sam- 
ple are available to the public for download from 
http : //www. astro . caltech.edu/ ovroblazars 



2. RADIO VARIABILITY 

To provide a robust measure of the amplitude of 
variability observed in a light curve, wc introduced 
the intrinsic modulation index, m = uq/Sq, where 
Co and Sq are the standard deviation and mean that 
would be measured for a light curve using a noiseless 
receiver with perfect sampling. In [B] we derive a like- 
lihood method for estimating m and its uncertainty 
from an observed light curve. Unlike the various mea- 
sures of variability amplitude commonly used in the 
literature, the intrinsic modulation index provides a 
rigorous measure of its uncertainty which is essential 
for statistical comparisons. 

We computed m for each CGRaBS source using two 
years of data in [6]. Using 3.2 years of data, we find 
that the intrinsic modulation index for each source 
is typically either consistent within uncertainties or 
is found to be larger with the longer data set. An 
examination of the sources with the most significant 
changes shows that the significant increases typically 
occur in sources that demonstrated no or little vari- 
ability during the first two years and subsequently ex- 
hibited a major flare or change in brightness. This 
suggests that the two-year data set was simply not 
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Figure 1: Likelihood distributions for the two-year mean 
intrinsic modulation index, mo, for the CGRaBS in 
1LAC (solid) and not in 1LAC (dashed). The gamma-ray 
loud 1LAC subset are 5.7% more variable on average, a 
> 6cr significant difference. 



long enough to capture the full range of variability 
behavior of the sources. Previous findings indicate 
that 4-6 year or longer timescales for flaring activity 
are common in blazars [e.g., [3], so further increases 
in the modulation indices are likely as the monitoring 
program continues. 



2.1. Population Studies 

In [6], we found strong (> 6a) evidence that 
the gamma-ray-loud subset of CGRaBS (i.e., those 
sources that appeared in the 1LAC catalog) was, 
on average, more variable at 15 GHz than was the 
gamma-ray-quiet subset (see Figure [I]). This trend 
persists in the longer data set with similar significance. 
This establishes the first rigorous demonstration of 
a connection between gamma-ray emission and radio 
variability amplitudes in blazars. 

In Figure [2] we compare the two-year variability 
amplitudes for CGRaBS sources identified as BL Lac 
objects with those identified as flat-spectrum radio 
quasars (FSRQs). The BL Lac objects are found to 
vary more strongly than the FSRQs with > 3er signif- 
icance. A similar result is found with the longer data 
set for CGRaBS sources. However, for the 1LAC sam- 
ple, using 3.2 years of data we find the BL Lac objects 
and FSRQs to be consistent with the same variability 
amplitude within ~ 1.5a (see Figure [3]). Further study 
of this effect in [7] suggests that the gamma-ray-loud 
FSRQs are much more radio variable than gamma- 
ray-quiet FSRQs whereas BL Lac objects display sim- 
ilar radio variability whether or not gamma-ray-loud. 
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Figure 2: Likelihood distributions for the two-year mean 
intrinsic modulation index, mo for the CGRaBS BL Lac 
object (solid) and FSRQ (dashed) subsamples. The 
BL Lac objects are 3.2% more variable on average, a 
> 3a significant difference. 
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Figure 3: Likelihood distributions for the two-year mean 
intrinsic modulation index, mo for the 1LAC BL Lac 
object (black) and FSRQ (red) subsamples. The two 
distributions are consistent with equal variability at the 
1.5cr level. 

2.2. Redshift Trend 

To examine the relationship between radio variabil- 
ity amplitude and redshift, we restrict our sample to 
the FSRQs. This is necessary because the lack of op- 
tical lines in BL Lac objects makes redshift determi- 
nation difficult, resulting in ~ 50% or less redshift 
completeness in our sample. This also prevents the 
difference in variability already found between BL Lac 
objects and FSRQs from contaminating our compari- 
son. In [B] , we found an apparent trend of decreasing 
variability amplitude with increasing redshift, shown 
in Figure [4] A comparison of bright (So > 400 mJy) 
FSRQs at z < 1 with those at z > 1 found the lower- 
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Figure 4: Average intrinsic modulation indices for 
CGRaBS FSRQs separated in redshift bins of width 
Az = 0.5. Error bars indicate the rms scatter within each 
bin. 



2.3. Cross-Correlation 

A major aim of the OVRO monitoring program is 
the detailed cross-correlation of radio light curves with 
gamma- ray light curves from the LAT. Peaks in the 
cross-correlation function could identify the relative 
locations of gamma-ray and radio emission within a 
blazar jet. However, establishing the physical signifi- 
cance of an apparent correlation is a major challenge. 
Typical methods for computing the significance of a 
cross-correlation peak do not account for correlation 
between points within each light curve and, as a result, 
seriously overestimate the significance of an appar- 
ent peak. We have developed a Monte Carlo method 
to evaluate the significance of apparent correlations 
more accurately. This method and early results are 
described in more detail in and in Max-Moerbeck 
et al. in these proceedings. A detailed study of cross- 
correlation between OVRO 15 GHz light curves and 
LAT gamma-ray light curves is forthcoming [5]. 
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Figure 5: Likelihood distributions for the intrinsic 
modulation index for low (z < 1, black) and high (z > 1, 
red) redshift subsets of the bright (So > 400 mjy) 1LAC 
FSRQ sample. The apparent difference is not 
statistically significant. 



3. KUPOL RECEIVER UPGRADE 

Blazars characteristically exhibit polarized radio 
emission which can reveal the configuration of the 
magnetic field in the synchrotron emission regions. As 
a result, monitoring radio polarization in addition to 
intensity can add valuable information and provide an 
additional channel to study the relationship between 
emission in the radio and other bands. To this end, we 
have designed and begun construction of a new Ku- 
band (15 GHz) receiver, "KuPol," with both linear 
polarization and total intensity measurement capabil- 
ities. Using a correlation receiver architecture and 
digital back-end, KuPol provides 16 MHz spectral res- 
olution over a total bandwidth from 12-18 GHz using 
the same beam-differenced optical design as the cur- 
rent program. Construction of KuPol is nearly com- 
plete and commissioning is planned for early 2012. 

4. CONCLUSIONS 



redshift FSRQs to be more variable with 3cr signifi- 
cance [5J. 

While this trend continues with the 3.2 year 
CGRaBS data, the significance of the difference is less 
with the additional data. We also find that FSRQs 
within the 1LAC sample do not show a significant dif- 
ference in variability between the high and low redshift 
subsamples (Figure [5]). In [7j, we find that this sig- 
nificance continues to decrease with additional data, 
suggesting that this observed difference is likely spu- 
rious. If real, however, a number of competing effects 
must be removed to determine whether the observed 
trend reflects any intrinsic variability trend with red- 
shift. 



We are actively and continuously monitoring a sam- 
ple of more than 1550 blazars, including all Fermi 
gamma-ray blazars north of 6 = —20° twice per 
week at 15 GHz. We have rigorously demonstrated 
a connection between gamma-ray emission and the 
strength of radio variability in our sample, an im- 
portant prerequisite for establishing the physical sig- 
nificance of correlations between radio and gamma- 
ray light curves. In addition, we have developed 
and begun applying a Monte Carlo method for ac- 
curately evaluating the statistical significance of ap- 
parent peaks in the discrete cross-correlation function 
between radio and gamma-ray light curves. In addi- 
tion to continuing our monitoring program, we will 
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soon deploy KuPol, adding an increased bandwidth, 
spectral resolution, and linear polarization monitoring 
to our radio monitoring program. These tools and the 
ever growing data sets will enable detailed study of the 
connection between radio and gamma-ray emission in 
blazars. 
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